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Abstract 

The  general  aim  of  this  project  is  to  get  an  improved  understanding  of  the  interac¬ 
tion  between  wall-generated  turbulence  and  compliant  surface  coatings  using  analysis 
and  direct  numerical  simulation  in  an  integrated  approach,  with  a  view  towards  the 
reduction  of  turbulent  sound  production  and  turbulent  drag. 

For  this  purpose,  in  a  first  step  that  is  targeted  at  identifying  interesting  domains 
in  the  space  of  parameters  describing  properties  of  a  compliant  wall  coating,  we  are 
developing  low-dimensional  models  based  on  Galerkin  projection  of  the  Navier-Stokes 
equations  onto  systems  of  eigenfunctions  obtained  via  Proper  Orthogonal  Decomposi¬ 
tion.  Because  of  the  relatively  small  effort  involved  in  simulating  and  analyzing  such 
models,  this  will  allow  us  to  scan  large  regions  of  parameter  space,  allowing  us  to  find 
regions  that  lead  to  a  reduction  of  turbulent  drag  and  turbulent  sound  production. 

Among  the  ultimate  goals  of  this  project  are  thus,  first,  to  obtain  a  fundamental 
understanding  of  flow-structure  interaction  phenomena  for  the  cEise  of  the  compliant- 
wall/turbulence  interaction,  and  second,  to  use  this  understanding  to  enhance  the  flight 
performance  of  air  vehicles  by  increasing  their  lift-to-drag  ratio. 
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1  statement  of  Objectives 

In  this  study,  research  is  conducted  into  the  issue  of  the  interaction  between  a  turbulent 
boundary  layer  and  an  adjacent  compliant  surface.  This  work  has  as  one  of  its  main  objec¬ 
tives  an  improved  imderstanding  of  the  dynamical  mechanisms  at  work  in  such  a  situation. 
Ultimately,  this  understanding  then  will  be  used  to  manipulate  the  turbulent  flow  such  that 
some  of  its  characteristics  axe  altered  in  a  desirable  way.  In  this  respect,  we  are  mostly 
interested  in  two  goals.  We  want  to  be  able  to  reduce  the  turbulent  drag  generated  by  such 
boundary  layers,  and  we  also  want  to  be  able  to  attenuate  the  noise  that  is  generated  by  the 
turbulent  flow. 

In  order  to  achieve  this,  we  will  try  to  And  regions  in  the  space  of  parameters  describing 
mechanical  properties  of  compliant  coatings  (stiffness,  damping)  within  which  the  interaction 
between  the  wall  and  the  turbulent  flow  is  such  that  the  above  goals  can  be  achieved. 

The  research  proposed  here  consists  of  two  main  components.  To  be  able  to  determine 
regions  in  the  parameter  space  of  the  wall  coating  that  are  of  interest,  we  will  construct 
low-dimensional  models  based  on  Galerkin  projections  onto  Karhunen-Loeve  eigenfunctions. 
These  models  can  only  approximately  describe  the  dynamics  of  the  turbulent  flow,  but 
because  of  their  low-dimensionality,  they  allow  us  to  look  at  Icirge  regions  of  the  parameter 
space.  Once  promising  combinations  of  parameters  are  determined,  we  will  then  use  accurate 
direct  numerical  simulations  to  assess  the  interaction  between  turbulence  and  the  compliant 
wall  in  more  detail. 


2  Status  of  Effort 

Our  work  on  developing  and  incorporating  models  for  the  interaction  between  the  turbulent 
flow  and  the  compliant  wall  is  proceeding  as  plainned.  As  a  reference  for  more  refined  models 
that  we  are  working  on  now,  we  have  developed  a  first  simplified  representation  of  the  wall 
as  a  three-degree-of-freedom  mass-spring  system  (see  below).  Results  from  a  study  of  the 
properties  of  that  model  are  encouraging,  and  are  presented  below.  A  refined  version  of 
this  model,  that  removes  most  of  the  limitations  of  the  previous  one,  is  currently  under 
investigation.  In  paxallel,  we  are  also  working  on  a  computer  code  for  the  direct  numerical 
simulation  of  turbulence/compli ant- wall  interaction.  That  code  will  probably  be  functional 
this  fall. 

3  Accomplishments 

3.1  Introduction 

In  the  following,  we  will  very  briefly  sketch  the  structure  of  our  models,  and  highlight  our 
most  important  results. 
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We  have  performed  an  analysis  of  the  boundary  conditions  (for  the  case  of  small  defor¬ 
mations  of  the  wall  as  in  the  study  of  Choi  et  al.  [3])  that  suggests  that  for  the  purposes  of 
our  low-dimensional  model,  neglect  of  the  nonlinear  terms  is  usually  quite  justifiable — the 
only  relatively  large  nonlinear  term  appears  in  a  term  for  the  streamwise  velocity  compo¬ 
nent,  which  makes  no  contribution  within  the  approximations  of  our  model.  The  boundary 
conditions  become: 

=  6 

U2  =  ^2  (1) 

^^3  =  6 

Note  that  if  the  velocity  field  is  expressed  as  a  linear  combination  of  eigenfunctions  Uj  = 
12n,k  exp(ifcz),  the  boundary  conditions  above  can  be  interpreted  as  a  restriction  on 

the  combinations  of  ^  and  xi  which  are  realizable. 

In  order  to  implement  these  boundary  conditions  in  our  model,  we  introduce  a  Fourier 
transform, 

Ui{x,t)=  (2) 

kiM 

where  we  can  also  write  u  using  our  POD  modes 

(3) 

n 

where 

am=  [  dy.  (4) 

Jo 

Neglecting  the  nonlinear  terms  in  the  waJl-normzd  boundary  condition,  we  obtain 

^2w  ~  ^2 

<^n<f>2w  —  ^2 

n 

0-x4>^2w  ^  ^2)  (3) 

where  it  is  assumed  that  the  first  eigenfunction  is  suflicient  to  represent  the  fluctuating 
velocity  field^ 

The  dynamical  equations  for  the  POD  amplitudes  (derived  from  the  Navier-Stokes  equa¬ 
tion  via  Galerkin  projection  onto  the  set  of  POD  modes)  become 

dm  +  •  •  •  =  4 — i>w~  +  •  •  • )  (6) 

P  o-i 

*  Note  that  the  “a”-sign  in  that  equation  has  to  be  taken  seriously:  At  this  level  of  truncation,  it  would 
in  fact  be  inconsistent  to  require  the  boundary  condition  to  be  met  using  just  one  POD  mode.  Observe  that 
Eq.  (5)  above  (with  a  “="  sign  instead  of  the  “w’’-sign),  once  ^2  a-nd  are  given,  uniquely  determines  a\. 
In  other  words,  the  dynamics  of  the  flow  would  be  directly  tied  to  the  dynamics  of  the  wall,  with  no  fluid 
mechanics  in  between.  Prom  a  dynamical  systems  point  of  view,  it  turns  out  that  we  need  at  least  one 
additional  degree  of  freedom  to  be  able  to  accommodate  both  the  weill  boundary  condition  and  the  Navier- 
Stokes  equation. 
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with  the  dots  indicating  the  terms  that  are  present  in  a  model  for  a  fixed  wall,  as  in  [1].  The 
pressure  at  the  wall  that  is  needed  above  can  be  obtained  by  using  the  instantaneous  Navier 
Stokes  equation,  which  finally  allows  us  to  write 


Pm  _ 

p  ikz 


(7) 


which  ultimately  gives 


Cm  +  •  •  • 


2- 


(8) 


It  may  be  a  bit  surprising  at  first  to  find  that  the  lateral  response  of  the  surface  apparently 
does  not  contribute.  The  reason  for  this  is  that — through  our  choice  of  considering  only  one 
POD  mode — the  boundary  conditions  at  the  wall  are  satisfied  only  in  an  approximate  fashion. 
This  is  also  the  reason  that  it  is  only  the  wall  normal  motion  that  affects  the  equations:  Once 
one  has  (implicitly)  decided  to  approximate  the  effect  of  the  wall  motion  on  the  boundary 
condition,  and  only  include  an  effect  through  the  pressure  term,  the  influence  of  tangential 
motions  is  gone  (since  in  the  pressure  integral  for  a  domain  with  a  horizontal  boundary 
only,  it  is  only  the  wall  normal  velocity  at  the  boimdaries  that  gives  a  contribution).  So 
it  is  the  one-mode  truncation^  that  causes  the  influence  of  tangential  motions  to  go  away, 
not  the  linearizations  that  we  have  introduced  above.  Note  that  this  approximation  is  quite 
consistent  with  the  overall  eiccuracy  of  this  type  of  model.  These  limitations  will  be  removed 
in  the  refined  model  described  below. 


3.2  Low- dimensional  Model  of  Flow  Over  a  Compliant  Wall 

3.2.1  Introduction 


We  have  investigated  a  model  for  a  turbulent  boundary  layer  interacting  with  a  deformable 
wall,  using  the  approximations  described  above.  The  flow  is  modeled  via  a  pair  of  complex¬ 
valued  ODEs  with  0(2)  symmetry,  and  the  dynamics  of  the  wall  is  described  using  a  simple 
damped  mass-spring  system.  The  complete  set  of  dynamical  equatiosn  consists  of  the  equa¬ 
tions  for  the  flow. 


dai 

dt 

dCL2 

dt 


^^(1) 

^lOi  -f  021^102  +  ai(eu|ail^  -f  ei2lG2p)  H - 


p.2^2  +  CnGj  -f  G2(62i|Q’iP  +  C22|a’2|^)) 


(9) 

(10) 


^In  conjunction  with  us  insisting  on  ai  being  determined  by  fluid  mechanics.  Note  that  if  we  had  gone 
the  alternative  route  of  interpreting  the  last  equation  of  Eq.  (5)  as  an  equation  for  ai,  then  the  dynamics 
of  the  flow  would  have  been  determined  entirely  by  the  wall  motions;  in  this  case,  we  would  have  to  decide 
which  component  of  the  wall  motion  to  pick,  because  prescribing  just  one  component  of  ^  would  already 
determine  ai . 
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and  a  second-order  equation  for  the  wall  motion, 

,(i) 

mi2  +  ci2  +  ki,=a,'^^.  (11) 

where  ^1“%^  =  0.2996aii/  (evaluated  numerically)  defines  the  pressure  at  the  wall  in 
Fourier  space. 

Figure  1  shows  a  comparison  of  the  behavior  of  this  model  for  the  cases  of  a  compliant 
and  a  fixed  wall,  respectively  (in  the  case  of  the  compliant  wall,  the  last  term  in  (9)  is 
ommitted,  together  with  (11).  One  can  see  that  for  the  small  value  of  c  (wall  damping)  that 


Hcieroclinic  Connection 


•2  1 


(a)  Fixed  wall. 


Wall  Nomol-Real  Subspace 


(b)  Compliant  wall.  Behavior  for  small  wall 
damping,  a  =  1.4,  i/  =  1.5,  u"  =  3,  c  =  1.0  • 
10“®,  fc  =  1. 


Figure  1:  Phase  portraits  of  low-dimensional  models. 

.was  chosen  in  the  figure,  the  lower  fixed  point  is  no  longer  attracting  in  all  directions  (it  is 
a  saddle  point).  Solutions  which  are  initially  attracted  to  this  point  oscillate  around  it  for 
an  extended  amount  of  time  before  diverging  to  infinity.  This  has  the  effect  of  reducing  the 
bursting  frequency  and  thus  potentially  reducing  turbulent  drag. 

3.3  Five  mode  model  with  wall-normal  compliance 

To  remove  the  limitations  of  the  model  described  above,  we  have  derived  a  refined  model, 
that  can  take  into  account  wall  motions  in  all  three  coordinates,  and  that  fully  satisfies 
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tlie  interface  conditions  between  the  moving  wall  and  the  flow.  To  do  this,  we  model  the 
compliant  surface  as  a  simple  mass-spring-damper  system  driven  by  the  fluid  stress  at  the 
wall.  In  the  streamwise  and  spanwise  directions,  the  deformations  of  the  wall  are  driven 
by  the  shear  stress  which  is  easily  accessible  in  the  models  as  a  linear  term  in  the  model 
coefficients  when  the  deformations  are  decomposed  into  Fourier  space: 

n 

=  fa*  Ui  =  E 

n 

(We  employ  capital  letters  to  describe  the  properties  of  the  comphant  surface  for  clarity.)  The 
wall-normal  deformations  of  the  comphant  surface  axe  driven  by  the  pressure  fluctuations  at 
the  boundary. 


(12) 

(13) 


Moll*  +  -^u;ll*  +  -f^iollfc  =  Pklwall  (14) 

Unfortunately,  the  pressure  at  the  boundary  is  not  as  easily  accessible  in  the  models  as  the 
shear  stress. 

The  pressure  at  the  wall  in  the  models  may  be  determined  by  examining  the  Galerkin 
projection  of  the  pressure  gradient  term  in  the  models. 

=  (pjc#:"’')  l“““  ^  (15) 

=  l-l  (16) 

The  eigenfunctions  are  defined  to  be  divergence  free  and  are  non-zero  only  close  to  the  wall, 
leaving  us  only  with  the  pressure  term  at  the  wall.  As  a  result,  the  model  equations  now 
have  an  additioned  term: 


^  ^  ^  ^  6.28a)6W"“')  a?'  +  ^ 

P  k',p,q 


r,k\p,q 


Our  previous  models  have  focused  exclusively  on  the  rigid-wall  case  with  (^2*^  =  0,  so  that  the 
pressure  term  disappears.  (If  the  eigenfunctions  are  defined  only  in  a  region  close  to  the  wall, 
an  additional  term  appears  which  was  incorporated  into  our  previous  models  as  a  forcing 
term.  The  magnitude  of  this  term  is  small  relative  to  the  others  and  will  be  neglected  here.) 
However,  to  incorporate  the  compliance  of  the  wall  and  satisfy  the  boimdary  conditions  in 
equation  1,  we  must  include  eigenfunctions  in  our  model  which  are  non-zero  at  the  wall. 

We  choose  to  base  our  low-dimensional  model  for  the  turbulent  boundary  layer  over  a 
compliant  surface  on  the  eigenfunctions  of  the  rigid-walled  boundary  layer.  We  will  intro¬ 
duce  additional  eigenfunctions  to  account  for  the  motion  of  the  compliant  surface  and  allow 
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the  boundary  conditions  to  be  satisfied.  In  the  absence  of  surface  compliance,  our  models 
will  revert  to  the  rigid-walled  case.  Since  we  do  not  have  an  experimental  or  computational 
database  on  the  turbulent  flow  over  a  compliant  surface  available,  we  must  derive  our  ad¬ 
ditional  eigenfunctions  in  an  ad  hoc  fashion.  We  choose  our  additional  eigenfunctions  as 
solutions  of  the  Stokes  equation  with  periodic  motion  of  the  wall: 

(I  ■ 

Uk{y  =  0,t)  =  cos{pt)  (19) 

where  beta  corresponds  to  the  natural  frequency  of  the  compliant  surface.  The  Stokes  equa¬ 
tion  for  the  streamfunction  if)  is  employed  to  determine  the  additional  eigenfunction  for  the 
wcdl-normal  motion  of  the  compliant  surface.  In  this  manner,  we  generate  three  additional 
eigenfunctions  —  one  for  each  direction  of  surface  motion  —  and  then  orthogonalize  them 
with  respect  to  each  other  and  the  rigid- walled  eigenfunctions  and  normalize  them. 

We  now  have  a  set  of  evolution  equations  for  the  coefficients  of  the  eigenfunctions  as  well 
as  evolution  equations  for  the  motion  of  the  wall  and  the  simplified  boundary  conditions. 
We  choose  to  determine  the  coefficients  for  the  additional  eigenfunctions  from  the  boundary 
conditions  at  the  wall.  If  only  one  or  at  most  two  eigenfunctions  has  a  particular  component 
which  is  non-zero  at  the  wall,  these  coefficients  may  be  found  easily.  The  pressure  may  then 
be  determined  from  the  evolution  equation  for  the  additional  eigenfunction  representing  the 
wall-normal  motion  of  the  wall 


Pk  lv=o  » 


-E 


(p) 

k 


(20) 


The  cubic  terms  disappear  because  this  additional  eigenfunction  is  defined  to  have  no  stream- 
wise  component,  and  there  is  no  streamwise  variation  in  our  model.  When  this  expression 
is  substituted  into  the  equation  for  the  wall-normal  motion  of  the  wall,  the  time  derivative 
term  merges  into  the  ^2*  term  resulting  in  an  new  effective  mass  of  the  wall  M^.  (In  fact, 
^2*  =  since  is  the  only  eigenfimction  with  a  non- zero  vertical  velocity  at  the 

wall.)  Normalizing  be  this  new  mass  results  in: 


^2*  d- 


fcp  0- 


(P) 

fc 


+ 


2-^  °‘k-k' 

k',P,q 


(21) 


Having  determined  the  coefficient  and  pressure  in  this  way,  we  have  a  complete  low-dimensional 
model  for  the  flow  over  a  compliant  surface. 


a(i)  - 


=E^’: 


kpO- 


ip) 

k 


+ 


k',p,g  r,k\p,q 


°‘rk’pq°‘k 


Re  (22) 


=  ElfcUaii  = 


(23) 
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Mifisfc  +  =  hk  \ivall  =  <ll”V3”2\ 


E6lM'”+E-=iVv),.<‘?’oK?<) 

P  fc'>Pi9 

(25) 


(26) 

(27) 

(28) 


We  are  in  the  process  of  using  this  low-dimensional  model  to  evaluate  the  potential  of 
compliant  surfaces  for  modification  of  the  dynamics  of  near- wall  turbulence  and  its  resulting 
effect  on  the  drag  at  the  wall.  In  addition  to  our  work  with  the  models,  we  have  begun  the  de¬ 
velopment  of  a  direct  simulation  code  for  simulating  turbulent  channel  flow  with  a  compliant 
boundary.  Because  of  the  complication  introduced  by  the  linearized  boundary  conditions,  we 
have  chosen  to  implement  the  compliant  boundary  using  the  immersed  boundary  technique 
which  was  developed  by  a  previous  member  of  our  group  [6]. 
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Division  of  Fluid  Dynamics,  Philadelphia,  PA,  November  1998.  Contributed. 

•  J.  Reuter  and  D.  Rempfer.  Direct  Niunerical  Simulation  of  Spatial  Pipe  Flow  Transi¬ 
tion  Using  a  Hybrid  Spectral/Finite-Difference  Scheme.  51st  Annual  Meeting  of  The 
American  Physical  Society’s  Division  of  Fluid  Dynamics,  Philadelphia,  PA,  November 
1998.  Contributed. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  Rice 
University,  Houston,  TX,  December  1998.  Invited. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near- wall  turbulent  flow.  United 
Technologies  Research  Center,  East  Hartford,  CT,  January  1999.  Invited. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  Applied 
Mathematics  Seminar,  University  of  Wisconsin  (Madison),  February  1999.  Invited. 

•  P.  N.  Blossey  and  J.  L.  Lumley.  Active  control  of  near-wall  turbulent  flow.  Dynamic 
Systems  and  Controls  Seminar,  University  of  California,  San  Diego,  March  1999.  In¬ 
vited. 

•  J.  L.  Lumley  and  P.  N.  Blossey.  Control  of  intermittency  in  near-wall  turbulent  flow, 
Isaac  Newton  Mathematics  Institute,  Cambridge,  UK,  June  1999.  Invited. 
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6.2  Technology  Transitions  or  Transfer 

None. 

6.3  New  Discoveries,  Inventions,  or  Patent  Disclosures 

None. 

6.4  Honors  /  Awsirds 

John  Lumley  is  a  Fellow  of  the  American  Academy  of  Arts  &  Sciences,  the  American  Physical 
Society,  the  American  Academy  of  Mechanics,  the  American  Institute  of  Aeronautics  and 
Astronautics,  and  was  a  Guggenheim  Fellow  from  1973-74. 

Dietmar  Rempfer  was  supported  by  a  Heisenberg  Grant  from  the  German  National  Sci¬ 
ence  Foundation  (Deutsche  Forschimgsgemeinschaft). 

Sidney  Leibovich  is  a  Fellow  of  the  American  Academy  of  Arts  &  Sciences  and  the 
American  Physical  Society. 
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